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Abstract 

Dynamical fermion mass generation is studied in the three-dimensional 
Thirring model reformulated as a gauge theory by introducing hidden local sym- 
metry. The analysis by use of Schwinger-Dyson equation is shown to exhibit 
a critical behavior as the number N of four-component fermions approaches 

2 



43 ■ ^cr = 128/3vH. 

> 

X 

When we consider the models of massless self-interacting fermions, a central issue 
is the dynamical fermion mass generation. If it is indeed the case, intriguing detailed 
questions are whether the model of interest is shown to exhibit a critical behavior 
as the number N of fermions approaches N cr and, in three dimensions, whether the 
pattern of dynamically generated masses preserves the parity or not. 

In this report we address ourselves to the problem of dynamical symmetry breaking 
in the Thirring model of N massless four-component fermions. (For the notations and 
the detailed calculations, see Ref.|l|.) 

G 

^Thi = Tpaij^d^lpa - — ^ ^Pal^a tpbl^b, (1) 
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where ijj a is a four-component Dirac fermion and the indices a, b run over N fermion 
species. Since the form of four-fermion interaction term is a contact term between vec- 
tor currents, a well-known technique to facilitate 1/iV-expansion is to rewrite Eq.([l|) 
by introducing an auxiliary vector field A 11 such as 

£ aux = E - 4=4*M. + ^V* M - ( 2 ) 



Gauge-noninvariant as the above Lagrangian is, however one may be tempted to 
regard as a gauge field. A systematic way to construct the £7(1) gauge theory, but 
is gauge equivalent to Eq.(|2|) is to elicit the fictitious Goldstone degree (or equivalently 
the Stiickelberg field) based on the principle of hidden local symmetry |J, J7], [j], 

£hls = E ^YD^ a + ^(A I1 - VNd^) 2 . (3) 

It is obvious that Eq.(^]) possesses a U{1) gauge symmetry and the gauge-fixed (uni- 
tary gauge) form of it exactly coincides with Eq.fQ), so does the original Thirring 
model in Eq. QT]) . 

Now that we find a gauge-invariant formulation of the Thirring model, we have the 
privilege to choose the gauge appropriate for our particular purpose. Here, instead 
of the unitary gauge notorious for loop calculations, let us consider the nonlocal 
gauge at the Lagrangian level 

£ GF = -\(d„A" + ^^)^gT ) (<W" + ViV^V), (4) 

where the gauge fixing parameter £ has the momentum- (derivative-) dependence. 
It is straightforward to prove the possession of the BRS symmetry in spite of the 
nonlocality of £ and thereby it guarantees the S-matrix unitarity. Another intriguing 
point is that, in the combined Lagrangian of Eqs. (^) and (||), the fictitious Nambu- 
Goldstone boson is completely decoupled independently of the specific form of 

In (2+1) dimensions, our Lagrangian in Eq.([3]) is invariant under the parity 

^ a (x) i-> ij/ a (x) = i 7 3 7 Va(x), A^x) ^ A'^x) = {-l) 5 ^A^{x). (5) 
and the so-called global "crural" transformation 

i> a ^ ^ = (exp(^ iQ y ® T a )i;) a , (6) 



2 



where S° = J, S 1 = — 27 s , X 2 = 7 5 , X 3 = — 7 5 7 3 and T a denote the generators 
of U(N). The question we shall address from now on is "which symmetry is broken 
dynamically?" In concern with the parity, first issue is whether one can take the 
regularization to keep both the U(l) gauge symmetry and the parity or not. Since 
our gauge action in tree level has the parity-conserving mass and the number of two 
component Dirac fermion species is even, the parity need not be violated by ap- 
propriate regulator. For example, the introduction of parity-conserving Pauli-Villars 
regulator leads to the parity-invariant effective action for the gauge field as have done 
in (2+l)D quantum electrodynamics(QED 3 ) [|IU|| . Another question is whether the 
pattern of dynamically-generated fermion mass involves the parity violating mass 
(— m^ a 7 5 7 3 '?/' a ) or not. Though at this stage we do not yet know whether the dy- 
namical symmetry breaking really occurs or not, in this gauge- invariant formulation 
of the Thirring model, such symmetry breaking pattern is proven to be energetically 



unfavorable by using the exact argument in Ref. II . Namely, since the tree-level 



gauge action corresponding to Eq.(|3|) is real and positive semi-definite in Euclidean 
space, energetically favorable is a parity conserving configuration consisting of half 
the 2-component fermions acquiring equal positive masses and the other half equal 
negative masses. 

According to the above arguments, the pattern of symmetry breaking we shall 
consider is not the parity but the chiral symmetry, i.e., of which the breaking is 
U(2N) — > U(N) x U(N). Thus we investigate the dynamical mass of the type mijjijj 
in the Schwinger-Dyson equation, giving 

(A(-p 2 )-lU-B(-p 2 ) = -- f-^-y A (-^ + B(-g 2 ) iD^(p-a) 
11 P) )P [ P) N J i{2^ 1 ^A\-q 2 )q 2 -B 2 {-q 2 ) u{P ' q) [P Qh 

(7) 

where the full fermion propagator is written as S(p) = i[A(—p 2 )^> — B(— p 2 )] -1 , and 
T u (p, q) and D^ijp — q) denote the full vertex function and the full gauge boson 
propagator, respectively. Task is, by employing some appropriate approximations, 
to reduce Eq.(^) to the tractable integral equation for the mass function M(—p 2 ) = 
B{— P 2 )/A(— p 2 ). First, we here adopt the 1/N expansion for T u (p, q) and D^ u {p — q) 
under a nonlocal gauge, i.e., they are the bare vertex and the one-loop vacuum 
polarization of massless fermion loop at the 1/N leading order. Then the Schwinger- 
Dyson equation (0) becomes the coupled integral equations for A(— p 2 ) and B(—p 2 ). 
They support a trivial solution A(— p 2 ) = 1 and B(—p 2 ) = at the 1/N leading order, 
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however, as was realized in QED3, we expect to find a nonperturbative nontrivial 
solution by examining them for finite N. 

A way is, by use of the freedom of gauge choice, to require A(— p 2 ) = 1 in a 
Schwinger-Dyson equation for A(—p 2 ). Then this gauge fulfills the consistency be- 
tween the bare vertex approximation and the Ward-Takahashi identity for the hidden 
local U(l) symmetry (or the current conservation), i.e., A(0) = 1. The specific form 
of the gauge is determined by a Schwinger-Dyson integral equation for A{— p 2 ), and 
it reduces the coupled Schwinger-Dyson equations into a single equation for B(—p 2 ) 
which turns out to be a mass function, i.e., M{— p 2 ) = B(—p 2 ): 



B{p 2 ) = ^f A d(q D ~ 2 )K(p, q] G)- 
J\ Jo ( 



■ q 2 + B 2 ( q 2 y 

where A is ultraviolet cutoff and the kernel K(p,q; G) is given by 

K(P, q\ G) = ^U^) [ dd sin D - 2 9 d{k 2 ) [D - V (k% (9) 
with k 2 = p 2 + q 2 — 2pqcos8. The gauge fixing parameter £ is a function of k 2 such 



as 



v(k 2 ) = a f^2 D '\t k2 (10) 



(D 



i(k 2 )G- i + k 2 

Gk D ~ 2 



\ / D 



D . Gk 1 
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where C^ 1 = T (2 — B (^-, ^ an< l 2-^1(0, b, c; z) the hypergeometric func- 

tion. Note that the kernel K(p, q; G) is positive-definite for positive arguments p, q, 
G and is symmetric under the exchange of p and q. 

Now let us show the existence of the nontrivial solution for the Schwinger-Dyson 
equation (jg) when 2 < D < 4. Since we have in mind the continuous phase transition, 
the solution of our interest is the nontrivial solution which starts to exist without gap 
in the vicinity of the phase transition point. Such a bifurcation point is identified by 
the existence of an infinitesimal solution 5B(p 2 ) around the trivial solution B(p 2 ) = 0. 
In terms of dimensionless variables (p = Ax 1 ^ ^ 2 ^, 5B(p 2 ) = AE(x), g = G/A 2 ~ D ), 
Eq.(||) is reduced to a linearized integral equation 

£(*) = I f 1 dyK(xW- 2 \yW- 2 h g)E(y), (12) 

iV Jam, 
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where o~ m = (m/A) D ~ 2 (0 < o m < 1) is the rescaled infrared cutoff and in fact m is 
nothing but the dynamically generated mass by the normalization m = 5B(m 2 ). We 
can rigorously prove that, if N is equal to the maximal eigenvalue of the kernel (0), 
then there exists a nontrivial solution E(a;). Hence, for a given a m , each line N(g, o~ m ) 
on (N,g) plane depicts a line of equal dynamically-generated mass m = Aa^ D ~ 2 \ 
Therefore, the critical line is defined by N cr (g) = \im Um _ N(g,o~ m ) which separates 
the broken phase from the symmetric phase. It is difficult to obtain the explicit form 
of the critical line N cr -(g) for arbitrary g, however we can get it in the limit of infinite 
four-fermion coupling constant, g — > oo. In this limit, the bifurcation equation ( |1"2"D 
in (2+l)D is rewritten into a differential equation 



d , 
dx 



,dY,(x) 



dx 
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(13) 



37T 2 iV 

plus the infrared boundary condition E'(cr m ) = and the ultraviolet one [xS'(x) + 
S(x) ] x= i = 0. When N > N cr = 128/37T 2 , there is no nontrivial solution of Eq.([L3|) 
satisfying the boundary conditions, while for iV < iV cr the following bifurcation solu- 
tions exist: 

/ x 



Or, 



sin(|5) W r , 



sm< 



{f 
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x 
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(14) 



where to = \j N cr /N — 1, 5 = 2lu 1 arctan to and a m is given by the ultraviolet bound- 
ary condition: 



V |T__ 



+ 25 



m 
A 



e 25 exp 



n 



1,2, 



(15) 
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The solution with n = 1 is the nodeless (ground state) solution whose scaling behavior 
is read from Eq.(|T5D: 

(16) 



y/N a /N-l- 

128/37T 2 is the same as the one in QED3 with 



The critical four-fermion number iV cr 
the nonlocal gauge. 

It is turn to comment briefly on the dynamically generated mass of the gauge 
boson and the dual transformation. The vector (gauge) boson is merely an auxiliary 
field at the tree level, however it turns out to be propagating by obtaining the kinetic 
term through fermion loop effect when the fermion acquires the dynamical mass. In 
(2+1) dimensions the pole mass My of the dynamical gauge boson is given by 



1 

2^ 



4m 2 + My 
2Mv 



V j. -1 My 



tan 



2m 



— m 



G~ 



(17) 
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which always satisfies a condition that My < 2m. Once the Thirring model is un- 
derstood as a gauge theory (|3]), its effective theory can be expressed in terms of dual 
antisymmetric-tensor field H lll ... /lD 2 of rank D — 2, which is actually a composite of 
fermions e^ 1 "'^ D d^ 2 H fJ , 3 ... flI3 = — TpT^/lf S a V'aT w V'a- This relation implies the current 
conservation at the quantum level, which is a crucial property for avoiding peculiarity 



in the lattice study of the Thirring model jTS). In (1+1) dimensions the above relation 



is nothing but the one of the bosonization, however extension to higher-dimensional 



case by integrating out the fermions leads to nonlocal bosonic Lagrangian flT2"fl . In 
(2+1) dimensions the dual gauge field shares exactly the same pole structure with 
the gauge field irrespectively of the phase. 

Analytic studies in the scheme of the Schwinger- Dyson equations with an 1/N ex- 
pansion yield the existence of the phase transition line on the plane of couplings (N, g) 
and predict the critical number of the four component fermion being N cr = 128/37T 2 
in the g — ► oo limit. In order to check these intriguing questions, an appropriate 
systematic method is lattice simulation. Here we present the status of the simulation 
of lattice gauge theory ]14[. The continuum Lagrangian is a U{1) gauge theory, and 
then the discretized version of Eq.(|3]) leads to a lattice gauge theory expressed in 
terms of a link variable e ie ^ and staggered fermions: 

where (3 = 1/g and 

M x , v = ^T,V^W W S y , x+ ^ - e-v^-ASy^) + m5 x , y . (19) 



Note that one can not discriminate whether this lattice Lagrangian flTB| ) at the tree 
level lies in the compact formulation or in the noncompact one since there is no gauge 
kinetic term at the tree level. This issue is rather subtle because each continuum 
limit may arrives at different continuum theory as has been done in QED 3 . 

So far, simulations for N = 2,4 and 6 systems on a 8 3 lattice volume with bare 
fermion mass m = 0.05, 0.025 and 0.0125 are completed by use of hybrid Monte Carlo 
algorithm and preliminary data on a 16 3 lattice for N = 2 and 6 are prepared for finite 
volume effect study |TJ| . Numerical simulation results on a 8 3 lattice are summarized 
in Fig. 1 where the plot points are those obtained by m — > extrapolation. N = 2 
shows smooth change of the chiral condensate and a long tail into weak coupling 
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regime, N = A shows a complex behavior of the condensate near the critical point, 
and N = 6 shows an abrupt change of the condensate which can be interpreted as a 
mean field type phase transition. The result from a coarse-grained lattice alludes to 
the existence of the critical fermion number N CI (2 < N CI < 6) for which the chiral 
behavior changes its character. It is consistent with the result in the continuum 
Thirring model reformulated as a gauge theory. Preliminary as the result is, the 
above suggests possible existence of the critical N CT in the lattice gauge formulation 
of the Thirring model. Data on a 16 3 lattice tells us that, though the dependence of 
the chiral condensate on both the size of the lattices and the extrapolation method is 
considerable quantitatively (particularly in the very weak coupling region of N = 2 
case), the overall qualitative nature of the chiral condensate does not seem to be 
changed. It is interesting to wait further lattice result of N = 6 case based on the 
auxiliary vecotr field [13j for the comparison with the above . 
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Figure 1: The chiral condensate, (ipip), vs. (3 for N = 2 (O points), iV 
points), and N = 6 (□ points). 
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In this report, we have reviewed how to reformulate the Thirring model as a gauge 
theory both on continuum and lattice, and, in such context, the dynamical symmetry 
breaking for iV four component Dirac fermions in (2+l)D. The Schwinger-Dyson 
approach in cooperation with 1/N expansion reaches a rigorous proof of the existence 
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of the chiral phase transition at a certain number of N and g under the light shed by 
the hidden local symmetry. The application of Vafa-Witten theorem forces the parity 
preservation in the dynamically-generated fermion mass spectra and the argument of 
parity-violating anomaly prevents the parity violation in gauge sector. In the infinite 
coupling limit (g — > oo), the Schwinger- Dyson equation was solved analytically and 
yields the critical fermion number N cr = 128/37T 2 in perfect agreement with N cr in 
QED3. More systematic analysis of the Schwinger- Dyson equations is in progess ||15|| . 
According to the present lattice simulations on 8 3 and 16 3 lattices, the lattice gauge 
formulation seems to support the existence of N cr (2 < N cr < 6), which should be 



tested again by further lattice simulation |T6 
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